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Abstract. We present a theoretical and experimental study of a quasi-one- 
dimensional zigzag antiferromagnet (N2H5)CuCl3, which can be viewed as weakly 
coupled Heisenberg chains with a frustrated interaction. We first discuss generic 
features of the magnetic properties of the zigzag spin chain between the nearly 
single chain case and the nearly double chain case, on the basis of the finite 
temperature density-matrix renormalization group (DMRG) calculations. We 
next show the experimental results for the magnetic susceptibility and the high- 
field magnetization of a single crystal of (N2H5)CuCl3 above the Neel temperature 
Tn = 1.55K. By comparing the experimental data with the DMRG results 
carefully, we finally obtain the ratio of the nearest and next-nearest exchange 
couplings as J1/J2 = 0.25 with J2/A;b = 16. 3K. We also investigate the three- 
dimensional (3D) coupling J' effect by using the mean-field theory combined 
with the DMRG calculations. The estimated value J' ~ O.O4J2 supports good 
one-dimensionality of (N2H5)CuCl3 above Tj^. 
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Figure 1. The zigzag spin chain. The dashed Hnes and sohd lines denote 
the nearest neighbor coupHng Ji and the next-nearest neighbor coupHng J2 
respectively. 



1. Introduction 

Frustration effects in low-dimensional antiferromagnetic (AF) quantum spin systems 
liave attracted much attention because of their peculiar behavior in low-energy 
physics. The quantum effect and the geometrical frustration induce strong fluctuation 
cooperatively, which often gives rise to various non-magnetic ground states and 
quantum phase transitions El E El ■ 

A typical example of such frustrated low-dimensional magnets is the = 1/2 
antiferromagnetic zigzag chain (see figure^. The Hamiltonian of the zigzag spin 
chain is given by 

^ -gtiBHj2s:, (1) 



Ti, — JiSi ■ S'i+i + J2Si ■ Si- 



where S is the 5*= 1/2 spin operator, g is the Lande's g factor and is the Bohr 
magneton. Ji and J2 denote the nearest neighbor and the next-nearest neighbor 
coupling constants respectively, and H is the applied magnetic field. Intensive 
theoretical studies have shown that the zigzag chain has a rich phase structure at the 
zero magnetic field. When J2 = 0, the system is equivalent to the Heisenberg chain 
having the gapless excitation [B]. Increasing J2, a quantum phase transition from the 
critical spin liquid phase to the gapped dimer phase |21[7| occurs at J2/J1 = 0.2411 jS]. 
This gapped dimer phase is expected to extend up to the limit J2/J1 = 00 which 
corresponds to the two decoupled Heisenberg chains. Here, it should be remarked that 
the dimer gap A for Ji ^ J2 is exponentially small: A ~ exp(— const x J2/J1) HI- 

The theoretical studies mentioned above have stimulated experimental studies of 
the zigzag spin chain. In fact, the magnetic properties of SrCu02 JOl and Cu[2-(2- 
aminomethyl)pyridine] 1111 were investigated recently. For the latter case, the nearly 
gapless behavior of the magnetic susceptibility was observed for a powder sample, 
and then the ratio of the couplings of Cu[2-(2-aminomethyl)pyridine] is estimated as 
J2/J1 ~ 0.2 |11| . In this paper, we focus on another type of the zigzag compound 
(N2H5)CuCl3 which is interestingly expected to have the coupling constant 

Ji ^ J2 CS|- Since (N2H5)CuCl3 can be synthesized as a single crystal, we can 
precisely compare its magnetic properties with numerical calculations. 

An important point on the compound is that its magnetic property is very similar 
to the one in the single Heisenberg chain limit, although it has the double chain 
structure. Thus how we can distinguish the magnetic behavior of the zigzag chains 
between Ji ^ J2 and Ji <C J2 is an important problem from both theoretical and 
experimental views. In what follows, we study the magnetic properties of a single 
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crystal sample of (N2H5)CuCl3, on the basis of the precise comparison between the 
numerical results obtained by the finite temperature density matrix renormalization 
group (DMRG) method (131 ^1 and the experimental measurements of the magnetic 
susceptibility and the high-field magnetization. In addition, we remark the effect of 
the three-dimensional (3D) couphng( J') on (N2H5)CuCl3; we apply the mean-field 
theory combined with the finite temperature DMRG to verify the one dimensionality 
of (NzHsjCuCla [131111. 

In section |21 we describe the finite temperature DMRG method and the 
fundamental aspects of the single crystal of (N2H5)CuCl3. In section O we show 
the DMRG results for the susceptibility and the magnetization curve, and discuss 
generic features of them. In section^ we show the experimental data of the magnetic 
susceptibility and the high-field magnetization parallel to the chain direction. Then 
the exchange coupling J2/A:b = 16. 3K and the ratio J1/J2 = 0.25 are obtained by the 
precise comparison of the experimental data with the DMRG results. We also discuss 
the 3D effect and estimate the coupling J' ~ O.O4J2. In section O the conclusion is 
summarized. 

2. Computational and experimental details 

2.1. finite temperature DMRG 

For the reliable analysis of the experimental results, systematic numerical data of the 
magnetization curve and the susceptibility are required in a wide parameter region. 
However, the quantitative calculation for the zigzag chain at finite temperatures has 
not been made yet especially for the nearly two chain cases. In this study we employ 
the finite temperature DMRG method ^J^], which is free from the negative sign 
problem and thus enables us to calculate the quantities sufficiently down to low 
temperatures with high accuracy. We note that the method is actually applied 
to some frustrated spin ladder systems successfully ^3 El 1201 . Following the 
implementation procedure described in reference I2()j . we calculate the magnetization 
and the susceptibility with the maximum number of the retained bases m = 240. The 
lowest temperature we can reach is fceT/ J = 0.07, where J is the larger one of the two 
exchange couplings Ji and J2. We have confirmed that the computed data converge 
with respect to m and the Trotter number. 

2.2. experiment 

We next summarize the synthesis and the crystal structure of (N2ll5)CuCl3. 
(N2H5)CuCl3 samples were synthesized according to the method described in 
reference |T21. After several syntheses, we obtained small single crystal samples with 
a typical size of 1.5 x 2 x 3mm^. Chemical analysis shows good agreement between 
the observed and calculated ratios of the elements. This compound crystallizes in 
the orthorhombic system (space Group: Pnma) J2|- The lattice constants at room 
temperature are a = 14.439(2)A, b = 5.705(1)1, and c = 6.859(1)1. The chain 
structure of this compound is shown in figure (21 This compound has the ladder 
structure in which Cu and CI align alternately along the leg, forming a copper zigzag 
chain. Copper 3d hole orbitals {3d^2_y2) are situated perpendicularly to the leg, 
resulting in small exchange interactions along the leg. 

Magnetic susceptibilities were measured with a SQUID magnetometer (Quantum 
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Figure 2. Chain structure of {N2H5)CuCl3. 

Design's MPMS-XL7) installed at KYOKUGEN in Osaka University. High-field 
magnetization measurements up to SOT were carried out with a pulse magnet at 
the same place. Besides, X-band ESR (Varian E109 spectrometer) and specific heat 
measurements were carried out at the same place. 

3. DMRG results 

In this section, we present the DMRG results for the susceptibility and the 
magnetization curve, and then discuss the fi'ustration dependence of such quantities 
for the nearly single chain and nearly decoupled chains cases. 

Figure |21 shows the normahzed susceptibilities x* = x/Xmax, where x is the 
magnetic susceptibility per site and Xmax is its maximum value. The normalized 
temperature T* = T/T^^i^-^ is also introduced, where Tmax is defined as the temperature 
at which x indicates Xmax- The values of Tmax and Xmax are summarized in table 1. 

Table 1. The values of Tmax and Xmax extracted from the DMRG results with 
the normalization Ji = 1 for Ji > J2 and J2 = 1 for Ji < J2. We also set 
/cb = SMS = 1 for simplicity. 





Ji > J2 






Ji < J2 




J2/J1 


^max 


Xmax 


J1/J2 


Tmax 


Xmax 


0.0 


0.640 


0.147 


0.0 


0.640 


0.147 


0.1 


0.588 


0.147 


0.1 


0.640 


0.143 


0.2 


0.516 


0.148 


0.2 


0.635 


0.139 


0.3 


0.419 


0.149 


0.3 


0.630 


0.135 


0.4 


0.299 


0.152 


0.4 


0.620 


0.131 


0.5 


0.279 


0.153 


0.5 


0.606 


0.128 



Let us first consider the nearly single chain case (Ji > J2). In figure 3(a), we can 




Figure 3. Finite temperature DMRG results of the susceptibility (a) for Ji > J2 
and (b) for Ji < J2. 



see the feature of the gapless Heisenberg chain for J2/J1 < 0.3; in a high temperature 
region x* increases with decreasing T*, and takes a round peak at T* ~ 1. Then 
X* starts to decrease, and finaUy converges to a finite value as T* ^ 0. Ahhough 
the spin fluid-dimer transition occurs at J2/J1 = 0.2411, the expected dimer gap 
is exponentially small for J2/J1 = 0.3 where the almost gapless behavior of x* 
appears 150]. As J2/J1 is increased to 0.4 and 0.5, the dimer-gapped phase starts to 
develop rapidly, and then the exponential decay of x* appears explicitly in the low 
temperature region. This rapid growth of the dimer gap also causes the considerable 
enhancement of x* in the high temperature side (T* > 1) for J2/J1 =0.4 and 0.5. 

We next consider the nearly decoupled chains case (Ji < J2), where we can also 
see that x* basically shows the feature of the gapless chain. A remarkable point in 
contrast to the nearly single chain case is that the gapless-like behavior is observed in 
the wide range of J1/J2 (see figure 3(b)) and thus the dependence of x* on J1/J2 is 
relatively weak. This wide-range gapless-like behavior is attributed to the extremely 
small spin gap: A ~ exp(— const x J2/J1), although the system is considered to be 
the dimer gapped phase for Ji ^ 0. For example, the spin gap is estimated as O.O3J2 
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even at J1/J2 — 0.5 and thus the effect of the gap is easily smeared out by the 
finite temperature effect. 

Here, we wish to point out an important feature of the susceptibihty when 
analyzing the experimental results. If the frustration is small, it may be difficult 
to discriminate between the behavior of the nearly single chain and that of the nearly 
decoupled chain only from the susceptibility result, because x* of the nearly single 
chain and that of the nearly decoupled chain are very similar to each other. However, 
an important point is that the frustration dependence of Tmax for Ji > J2 is relatively 
sensitive to J2/ Ji(see table 1), implying that the energy scale of the dominant coupling 
constants has a certain difference between Ji > J2 and Ji < J2, even if the shapes 
of the susceptibility are quite similar to each other. For example, the shape of x* for 
J2/J1 — 0.1 is almost coincident with that for J1/J2 = 0.3, but the values of T^ax for 
J2/J1 — 0.1 and J1/J2 = 0.3 can read as 0.588 and 0.630 respectively, which show 
about 7% deviation. 

We now proceed to the results of the magnetization curves. Figure 4 shows 
the magnetization curves both for Ji > J2 and for Ji < J2 at fcBT'/Ji(J2)=0.125, 
which can be directly compared with experimental results of (N2H5)CuCl3 measured 
at 2K as shown later. In the figure, we can basically see the characteristic property 
in the gapless Heisenberg chain as well; after the linear behavior with H is observed 
near the zero magnetic field {H — 0), the magnetization M increases rapidly in the 
high field region, and finally saturates near gfJ-sH/ Ji(J2) ~ 2.5. However, what we 
want to remark here is that the quantitative differences appear in the magnetization 
curve. For instance, in figure 4(a), we can see that the growth of the dimer gap for 
J2/J1 > 0.4 causes the enhancement of the magnetization in the middle-field region 
{g^sH/Ji ~ 1), which is consistent with the susceptibility result. Moreover the 
magnetization curve around the saturation field Hg is shifting to the high field side 
systematically as the frustration increases. 

In order to capture the feature of the magnetization curves in each region of 
Ji > J2 and Ji < J2, we here focus on the above mentioned systematic shift of the 
magnetization curve near the saturation field Hs. For the zigzag chain, the behavior 
of the saturation field Hg is easily obtained as 



For the nearly single chain, the saturation field Hg stays at 2Ji/g^B for J2/J1 < 1/4, 
and it starts to increase gradually as J2 / Ji increases beyond 1 /4 . On the other hand, 
for the nearly decoupled chain, Hg starts to increase linearly as soon as the inter-chain 
coupling is introduced. This dependence of Hg readily appears in the magnetization 
curve shown in figure^l the magnetization curves near Hg for J2/J1 < 0.3 remains at 
the "vicinity" of the pure Heisenberg one. In contrast, for the nearly double chain, 
we can see that the curves near Hg shift to the high field side linearly with respect 
to Ji/J2- As a result, the magnetization curves near Hg for Ji > J2 become slightly 
steeper than those for Ji < J2. 

As was seen for the susceptibility result, the dominant coupling determined from 
the susceptibility exhibits a slightly different value depending on the nearly chain or 
nearly decoupled chains cases. Then such difference of the dominant coupling affects 
the magnetization curve near the saturation field, which can be resolvable with the 




for (J2/J1 < 1/4) 



for (J2/J1 > 1/4) . 
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comparison between the accurate DMRG result and the experimental ones. Therefore 

the coupling constant of the zigzag chain can be expected to be determined with the 
cooperative use of the susceptibility and the magnetization curve. 




Figure 4. The magnetization curves obtained with the DMRG (a) for Ji > J2 
at knT/Ji = 0.125 and (b) for Ji < J2 at k^iT/J^ = 0.125. 



4. Experimental results 

4-1- ESR and specific heat 

We performed X-band ESR measurements at room temperature on a single crystal 
sample of (N2H5)CuCl3 in order to get the g'-values precisely. Accurate comparison 
between the experimental data and the numerically calculated ones rcqiiircs the precise 
g-values. Figure 5 shows the angular dependence of the g-value. We obtained the g- 
values of 2.285 and 2.060 for the external magnetic field parallel and perpendicular 
to the chain, respectively. In the following comparisons between experiments and 
calculations, we used these ^-values. 
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Figure 5. The angular dependence of the g-value of (N2H5)CuCl3. 
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Figure 6. The specific heat of a powder sample of (N2H5)CuCl3. 



We carried out specific heat measurements on a powder sample of (N2H5)CuCl3. 
Figure 6 shows the temperature dependence of the specific heat divided by 
temperature. Here, the data include the lattice contribution. We observed an 
anomaly at 1.55K (Tn) which is probably due to the antiferromagnetic long-range 
ordering. In order to avoid the 3D effect in this quasi-one-dimensional system, high 
field magnetization experiments were done at 2K above Tn. 

4-2. magnetic susceptibility 

Magnetic susceptibility Xbxp of the single crystal of (N2H5)CuCl3 parallel to the 
chain, which was used for the ESR measurement, is presented in figure [3 where the 
DMRG results Xdmrg are also shown for comparison. In the high temperature region 
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Xexp exhibits the Curie- Weiss-type T-dependence, and has a round peak around 
T = lOK. When the temperature is decreased further, the Xexp also decreases and 
finally converges to a finite value in the zero temperature limit. Since these features 
are characteristic in gapless ID antiferromagnets, we should take account of the two 
possibilities of the gapless zigzag chains, as mentioned in the introduction. We thus 
analyze Xexp with the DMRG results both for the single chain limit (Ji > J2) and 
the weakly coupled Heisenberg chain limit (Ji < J2). In fitting, we first estimate the 
dominant exchange coupling constant from the position of the broad peak, by using 
(7=2.285 obtained from the ESR, and then tune another exchange coupling constant 
so as to reproduce the shape of the susceptibility well. 

For the former case (Ji > J2), we find tha^t Xdmrg 

of J2/J1 = 0.1 with 

Ji/ku =17.4K agrees well with Xexp in the whole temperature region (see figured (a)). 
For the latter (Ji < J2), we can also see that the DMRG results for 0.2 < J1/J2 < 0.3 
well reproduce the shape of the experimental susceptibility as in figure[7|(b), where the 
best fitted value is obtained as J1/J2 = 0.25 with J2/kB = 16. 3K. Since the DMRG 
results seem to explain the whole shape of the susceptibility well for both cases, it is 
unfortunately difficult to confirm which case is realized in (N2H5)CuCl3 only from the 
susceptibility result, as discussed in the previous section. However, we note that the 
dominant couphng Ji/ks =17. 4K for Ji > J2 is about 6% bigger than J2/kB — 16. 3K 
for Ji < J2, which is consistent with the DMRG result. For the case oi H 1. chain, 
very similar results were obtained using the same parameter values of the exchange 
couplings as for H / / chain and (7=2.060 from ESR. 

Before proceedings to the next subsection we make a comment about a possibility 
of the weak ferromagnetic coupling. We have also computed the case of a negative 
J2{Ji) for Ji > J2(Ji < J2)- However, the shape of the susceptibihty is changed 
drastically by introducing the ferromagnetic coupling, and so that the couplings of 
(N2H5)CuCl3 are antiferromagnetic. 

4-3. magnetization curve 

Figure |HI shows the magnetization curve for the single crystal of (N2H5)CuCl3 
along the chain at T = 2. OK up to 30T. The temperature was set a little above 
Tn (1.55K) as already mentioned. Around the saturation field, there is a small 
hysteresis of magnetization between the ascending and descending processes, due 
to the magnetocaloric effect. In the figure, we can see that the measured curve of 
(N2H5)CuCl3 is qualitatively similar to the gapless chain, which is consistent with 
the susceptibility results. Using the parameters estimated in the previous subsection, 
we now discuss the magnetization processes of the DMRG and the experiments. The 
temperature used in the numerical calculations was determined close to that in the 
experiments. 

The results of comparisons are also shown in figure |S| where the DMRG curves 
both for Ji > J2 and for Ji < J2 are in good agreement with the experimental one 
in the low-field region. However, the magnetization curve in the high field region is 
reproduced well by the coupling of J1/J2 = 0.25, as can be seen in figure|S^b). Indeed 
the gradient of the calculated curve for Ji > J2 is slightly bigger than the experimental 
one around 15 < H < 25T, while the curve for Ji < J2 is well fitted in the all range of 
the magnetic field. The difference between two cases described in the previous section 
certainly appears in high-field region. Hence we have determined that the zigzag chain 
compound (N2H5)CuCl3 is a weakly coupled double chain system with the parameters 
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Figure 7. Magnetic susceptibility of {N2H5)CuCl3. The open circles denote 
experimental data Xexp and the lines show finite temperature DMRG results 
Xdmrg- 



J1/J2 = 0.25 and J2//sb = 16.3K. 
3D interaction 

In this subsection, we would like to discuss the effect of the 3D interaction on 
(N2H5)CuCl3. Since the specific heat shows the 3D ordering at Tn = 1.55K, some 
consideration for the 3D coupling may be required. To this end, we employ here the 
mean-field theory combined with the DMRG |17[ [TH] . We use the obtained values of 
the couphngs J1/J2 = 0.25 and J2/A;b = 16. 3K, and assume a 3D Heisenberg type 
interaction whose coupling constant is denoted as J'. According to the strong leg 
interaction J2, we assume the staggered magnetic field along the leg direction, and 
then determine the assumed magnetic field and the DMRG-calculated magnetic order 
self-consistently. In particular, the phase transition point within the mean field theory 
is given by 



1 



(3) 



where Xs{T) 
Calculating 



is the staggered susceptibility and z is the coordination number. 
J with the finite temperature DMRG, we can estimate J' from the 
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Figure 8. The magnetization curves of (N2H5)CuCl3. The comparison with the 
DMRG results (a) for J2/J1 = 0.1 and (b) for J1/J2 = 0.25. 



experimentally determined Neel temperature Tn through the equation Q- In figurc|3 
we plot the normalized staggered susceptibility = J2XS obtained with the finite 
temperature DMRG. In the figure, we find that the crossing point of xt a-nd the 
experimentally obtained Neel temperature is given by xt — 6-6 and = k^T^/ J2 ~ 
0.095, where the T-^ is the normalized transition temperature. Assuming 2=4, we 
obtain J' ^ 0.04J2('^ 0.6K), which is about 0.2Ji. Thus we can conclude that the 3D 
interaction effect is much smaller than the intra-chain couplings. Therefore the zigzag 
chain picture is valid for (N2H5)CuCl3 above the Neel temperature. 

5. Summary and discussion 

We have investigated magnetic properties of (N2H5)CuCl3, which can be regarded as 
a zigzag chain with the coupling Ji < J2 by the accurate comparison between the 
numerical calculations and experiments. 

First, we have calculated the susceptibility and the magnetization curve of 
the zigzag chain using the finite temperature DMRG method for the nearly single 
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Figure 9. Staggered susceptibility of the zigzag chain for J1/J2 = 0.25. The 
dotted Hne denotes the normaUzed transition temperature = k-BTfi/J2 = 
0.0951 ■ ■ ■ for (N2H5)CuCl3. The value of xt at k^T j J2 = gives the coupling 
zj' following equation ISl . 

chain case Ji > J2 and for the nearly double chain case Ji < J2 to clarify the 
differences between the two cases appearing in the parameter dependence of the 
physical quantities. From the experimental point of view, the quantitative difference 
appearing in the high-field region of the magnetization curve is found to help us to 
identify the magnetic nature of (N2H5)CuCl3. 

We have next presented the experimental results for the single crystal of 
(N2H5)CuCl3. The values parallel and perpendicular to the chain are determined as 
g — 2.285 and 2.060, respectively by the ESR experiment. The Neel temperature is also 
estimated as Tn = 1.55K from the specific heat of a powder sample. Above the Neel 
temperature, the observed susceptibility and the magnetization for the single crystal 
show gapless or nearly gapless behavior. An accurate comparison of the experimental 
data with the DMRG results shows that the DMRG results for J1/J2 = 0.25 with 
'/2/fcB = 16. 3K well reproduce the experimental results. The estimated couplings also 
show that the dimer gap of (N2H5)CuCl3 is O.OOIJ2 ^ 0.02K 9', which is much smaller 
than Tn. Hence the gapless behavior of (N2H5)CuCl3 observed above Tn is consistent 
with the obtained exchange couplings. 

Finally, we have discussed the effect of the 3D interaction, using the mean-field 
theory combined with the finite temperature DMRG. The estimated value of the 3D 
coupling J' (~ O.O4J2) is sufficiently small compared with the intra-chain couplings, 
which supports the zigzag chain picture of the compound. 

In the present paper, we have addressed the zigzag compound of (N2H5)CuCl3 
in terms of the bulk quantities: the magnetic susceptibility and the high-field 
magnetization, for which the rather precise analyses are required. However, the 
frustration effect often gives rise to the incommensurate behavior in the correlation 
function which is difficult to be captured by the bulk quantities. Thus it 

is also an interesting problem to investigate the compound by measurements of the 
fluctuation-associated quantities with the position resolvable methods, such as NMR 
or the inelastic neutron scattering. 
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